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After the first observation of the double giant dipole 
resonance (DGDR) in relativistic heavy ion coUisions 
(RHIC) l^g] the magnitude of its excitation cross sec- 
tion attracts a special attention. This is because the first 
data j2j indicated a strong enhancement of the DGDR ex- 
citation in ^^^Xe (by a factor of 2-4) as compared to the 
theoretical predictions available Several attempts 

have been made to understand the reason of this phe- 
nomena 1^-^. A few years later a similar experiment 
with the excitation of the DGDR in ^°*Pb brought an- 
other intriguing news [ pO| which were interpreted as a 
good agreement between experiment and theory if some 
corrections to the pure harmonic picture of the DGDR 
excitation were added |^-|^. Unless a good systematic 
study is achieved, the question on the excitation cross 
section of the DGDR in RHIC remains open. This stim- 
ulates theoretical studies on the different processes which 
might be responsible for an enhancement of the DGDR 
strength function. 

Recently it was argued that transitions between com- 
plex configurations, to which the GDR and DGDR door- 
way states are coupled, may be responsible for the en- 
hanced DGDR decay into the GDR states as compared 
to the GDR decay into the ground state However, 
it should be remembered that because of the available 
phase space two consequent 7-emissions from the DGDR 
is not the same as the inverse process of the DGDR exci- 
tation. The transitions between the GDR and the DGDR 
complex configurations were not taken into account in 
previous microscopic studies |]9|, p^ , |l3| . The role of these 
transitions for the DGDR excitation in RHIC will be con- 
sidered in the present paper. It will be concluded that 
their role is marginal in the process under consideration 
although a huge amount of the i?l-strength is hidden 
in the GDR^DGDR transition. This negative result en- 
sures that calculations, in which only transitions between 
collective components of the GDR and DGDR are taken 
into account and which are much easier to carry out, re- 
quire no further corrections. 

The main mechanism for the DGDR excitation in 
RHIC within a semiclassical approach is a two-step pro- 
cess g.s.^GDR^DGDR |l4|. Corrections to the second- 
order perturbation theory arising from coupled-channels 
were studied in Ref. j^]. Although for grazing impact 
parameters the coupled-channels calculation deviates on 
the 20% level from second-order perturbation theory, it 
makes only a small change to the total cross sections. 
That is, for not too small impact parameters the second- 
order perturbation theory works quite well. Indeed, this 
has been observed recently in the analysis of the experi- 
ment on DGDR excitation in lead projectiles impinging 
on different Z-targets |]l0| . 

Thus, we can use for the excitation probability of the 
DGDR 
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where the index i labels intermediate states belonging to 
the GDR, and 

^Ji(Mi)^j2(M2) first-order El exci- 

tation amplitude for the transition Ji(Afi) — > J2(-^^2) 
in a collision with impact parameter b. For each state, 
J and M denote the total angular momentum and the 
magnetic projection, respectively. 



The amplitude aj^j-^/j^-j 



+ J2(M2 



is given by 

"^,(M,)^J2(M2)(^' ^) = {JlM.MJ^M^) 

X <J2\\El\\Ji> fEl(^.){E, b) . 

It is a product of the reduced matrix element < 
J2||i?l||Ji > for the i?l-transition between the states 
Ji(Mi) and J2(-^2) which carries nuclear structure in- 
formation and the reaction function fEi(n){E^ b). The 
latter depends on the excitation energy, charge of the 
targ et, beam energy, and is calculated according to Ref. 
||l5| . Except for the dependence on the excitation energy, 
it does not carry any nuclear structure informati on. T he 
cross section for the DGDR is obtained from Eq. (0.1) by 
integration over impact parameters, starting from a min- 
imal value bmin to infinity. This minimal value is chosen 
according to Ref. 

In microscopic approaches the strength of the GDR 
is split among several one-phonon 1~ states (due to the 
Landau damping). The wave function |1~ > couples to 
complex configurations |1^ > yielding the GDR width. 
We use the index a for simple configurations and the 
index (3 for complex ones, respectively. Thus, the wave 
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where coefficients S^°^{a) and C!^°^{(3) can be obtained 
by diagonalizing the nucl ear model Hamiltonian on the 
set of wave functions (0.2). 

The total i?l-strength of the GDR excitation from the 
ground state, Boor(^I) = I < irll^l||0+.. > ?, re- 
mains practically the same as in the one-phonon RPA 
calculation because the direct excitation of complex con- 
figurations from the ground state is a few order of mag- 
nitude weaker as compared to excitation of one-phonon 
states. However these complex configurations play a fun- 
damental role for the width of the GDR. 

The wave function of the 2+ component of the DGDR 
states can be written in the similar fashion: 

a — {a.i X a2 } 
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In this equation we separated in the first term the door- 
way [1~ X 1~] DGDR configurations from other two- 
phonon configurations (second term) and complex con- 
figurations (the last term). The same equation as (0.3) 
is valid for the 0+ DGDR states (see Ref. Q for a dis- 
cussion on the role of the 0+, 1+, and 2+ DGDR states 
in the excitation process). 

The total i?l-transition strength between the GDR 
and DGDR, E/E, I < 2+(0+)/||i;i||l- > is much 
larger as compared to that for the GDR excitation, 
J2i \ < 1,"II-E1||0+^. > |2, from the ground state. This 
is because the former includes transitions not only be- 
tween doorway GDR and DGDR states but also between 
complex configurations as well. The enhancement factor 
should be the ratio between the density of doorway and 
complex configuration in the GDR energy region. But in 
the two-step excitation pro cess the sum over intermedi- 
ate GDR states in Eq. ( |D.l[ ) reduces the total transition 
strength for g.s. GDR DGDR to - 2 • |Bgdr(-B1)P 
(the factor 2 appears due to the bosonic character of 
the two phonons which also holds if Landau damping 
is taken into account). To prove this we substitute the 
wave functions of the GDR and DGDR states given by 
Eqs. (0^,3.3) in expression (p.l[). We obtain two terms. 
The first one corresponds to transitions between doorway 
GDR and DGDR states (after the GDR is excited from 
the ground state through its doorway component): 

^M.' = E E b) < i^ii£;iiio+,. > 

i a a' a. 
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and the second one accounts transitions between complex 
configurations in the wave functions of Eqs. (|D^J0|): 



proportional to 2 • |i?GDR(-£'l)P in second order perturba- 
tion theory. This argument was the reason for neglecting 
the term i?^^' in previous calculations of DGDR exci- 
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The second reduced matrix element in the above equa- 
tions is proportional to the reduced matrix element be- 
tween the ground state and the doorway one-phonon con- 
figuration p. 

For a given impact parameter 6, the function 
fEi(p){E, h) can be approximated by a constant value 
fsiifj.) [0 relevant values of the excitation en- 

ergies. Then the energy dependence can be taken out 
of summations and orthogonality relations between dif- 
ferent components of the GDR wave functions can be 
applied The orthogonality relations between the 

wave functions imply that J2i 5f °^(a)Cf °^(/3) = 0. This 
means that the term B^^i vanishes. The term A^^^i 
summed over projections and all final states yields a tran- 
sition probability to the DGDR, PoGDRiEf , b), which is 



tation in Rcfs. |y,|12|,y_3| where the coupling of doorway 
GDR and DGDR states to complex configurations was 
taken into account. 

In Fig. [l] we plot the value of XEiiE) = 
2tt f db bJ2f_i \fEi{^L){E, &)P as a function of energy cal- 
culated for the 208pb (640-A MeV) 4- ^ospb reaction. 
This value corresponds to CTgdr if Sgdr(£'1) = 1- The 
square in this figure indicates the location of the DGR 
in ^°*Pb. This figure demonstrates that the function 
Xei{E) changes by 60% in the GDR energy region. The 
role of this energy dependence for other effects has been 
considered in Refs. Taking into account that one- 

phonon 1^; configurations are fragmented over a few MeV 
|l3| , when a sufficiently large two-phono n b asis is in- 
cluded in the wave function given by Eq. (0.2), the role 



of the Bf^fj^i term in the excitation of the DGDR should 
be studied in more detail. 

To accomplish this task we have performed firstly a 
simplified calculation in which we used the boson type 
Hamiltonian: 
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where Qj, is the phonon creation operator and uoa is the 
energy of this one-phonon configuration; is the oper- 
ator for creation of a complex configuration with energy 
Zjfi and Up is the matrix element for the interaction be- 
tween these configurations. We have assumed that the 
energy difference between two neighboring one-phonon 
configurations is constant and equals to Aw. An equidis- 
tant spacing with the energy Aw was assumed for the 
complex configurations. We also have used a constant 
value U for the matrix elements of the interaction. The 
Sgdr(-£'1) value was distributed symmetrically over door- 
way one-phonon configurations. Thus, the free parame- 
ters of this model are: Aw, Aw, C/, the number of one- 
phonon and complex configurations, and the distribution 
of the i3GDR(^'l) value among the doorway states. The 
only condition we want to be satisfied is that the energy 
spectrum for the GDR photoexcitation is the same as the 
one known from the experiment. 

After all parameter s are fixed we diagonalize the model 
Hamilto nian of Eq. ( |0.6[) on the set of wave functions 
of Eq. (|0^) for the GDR and on the set of Eq. (joj) 
for the DGDR. The diagonalization procedure yields the 
information on eigen energies of the 1~ GDR states and 
on the coefficients Sf°^{a) and Cf°"(/3), respectively. 
One also obtains information on eigen energies of the 2jr 

or Oj^ DGDR states and the coefficients S°f'''°^{a) and 
(^DGDR^^/-j^ respectively. With this information we are 
able to study the role of the B^^i term in the excitation 
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which cTogdr I^) = o-dgdr(£^) ^ Jdb 
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of the DGDR in RHIC. 

The big number of free parameters allows an infinite 
number of suitable choices. In fact, not all of the param- 
eters are really independent. For example, the increase 
in the number of simple or complex configurations goes 
together with the decreasing of the value of U. This is 
necessary for a correct description of the GDR photoab- 
sorption cross section. This makes it possible to inves- 
tigate the role of the term in different conditions 
of weak and strong Landau damping and for different 
density of complex configurations. In our calculations 
we vary the number of collective doorway states from 
one to seven and the number of complex configurations 
from 50 to 500. The value of U then changes from about 
100 to 500 keV. The results of one of these calculations 
for the excitation of the 2+ component of the DGDR in 
2°spb (640-A MeV) + ^"Spb collisions are presented in 
Fig. H For a better visual appearance the results are 
averaged with a smearing parameter equal to 1 MeV. 
The dashed curve shows the results of a calculation in 
A 'J^^ I' ^bb\Af,f,,\^ and there- 

^ '^dgdr(-^) ~ 

J dbb \A^^> + B^^'l'' are represented by a solid curve. 

Our calculation within this simple model indicates 
that the role of the Bfj_fj_' term in the second order per- 
turbation theory is negligibly small, although the total 
B{E1) strength for transitions between complex GDR 
and DGDR configurations, considered separately, is more 
than two orders of magnitude larger than the ones be- 
tween doorway GDR and DGDR configurations. The 
value Act = (cToc^d^ - f^DGDR)/o'DGDR: where a^Ht^^ = 
/ <^oGTm^\E)dE , changes in these calculations from 1% 
to 2.5%. The results practically do not depend on the 
number of complex configurations accounted for. The 
maximum value of Act is achieved in a calculation with 
a single doorway GDR state (no Landau damping) . This 
is because the value of U is the larger in this case and the 
fragmentation of the doorway state is stronger. Thus, in 
such a situation, the energy dependence of the reaction 
amplitude modifies appreciably the orthogonality rela- 
tions. But in general the effect is marginal. 

We also performed a calculation with more realistic 
wave functions for the GDR and DGDR states taken from 
our previous studies |^] which were based on microscopic 
QPM (Quasiparticle Phonon Model) |l|]. These wave 
functions include 6 and 21 doorway states for the GDR 
and DGDR, respectively. The complex configurations are 
two-phonon states for the GDR and three-phonon states 
for the DGDR. The energies of the doorway states and 
complex configurations were obtained from RPA equa- 
tions and the matrix elements of the interaction, , were 
calculated on a microscopic footing without free param- 
eters by making use of the QPM Hamiltonian and inter- 
nal fermion structure of phonons. The value Act equals 
in this realistic calculation to 0.5%. This result is not 
surprising because realistic calculations with only two- 
phonon complex configurations, and a limited number of 



them, somewhat underestimate the GDR width which is 
crucial for the modification of the orthogonality relations. 

In conclusion, we investigated the role of transitions 
between complex GDR and DGDR configurations within 
second-order perturbation theory for the DGDR excita- 
tion in RHIC. We have proved that these transitions play 
a marginal role in the process under consideration and it 
is sufficient to take into account only transitions between 
the ground state and doorway GDR and DGDR config- 
urations. 

V.Yu. P. thanks for the hospitality at the Insti- 
tuto de Ffsica of the Universidade Federal do Rio de 
Janeiro, where this research has been performed, the 
CNPq for the financial support and Prof. P.F. Bor- 
tignon for fruitful discussions. This work was par- 
tially supported by the RFBR (grant 96-15-96729), by 
the Heisenberg-Landau program, by the FUJB/UFRJ 
and by the MCT/FINEP/CNPq(PRONEX) (contract 
41.96.0886.00) 



[2: 
[s: 
[4: 

[6; 

[7: 
[s: 
[9: 
[10: 

[11 

[12 

[13 

[14 
[15 
[16 



J.L. Ritman et al. (TAPS collaboration), Phys. Rev. Lett. 
70, 533 and 2659 (1993). 

R. Schmidt et al. (LAND collaboration), Phys. Rev. Lett. 

70, 1767 (1993). 

C.A. Bertulani, and V.G. Zelevinsky, Phys. Rev. Lett. 

71, 967 (1993); Nucl. Phys. A 568 931 (1994). 

V.Yu. Ponomarev, E. Vigezzi, P.F. Bortignon, R.A. 
Broglia, G. Colo, G. Lazzari, V.V. Voronov, and G. Baur, 
Phys. Rev. Lett. 72, 1168 (1994). 

C.A. Bertulani, L.F. Canto, M.S. Hussein, and A.F.R. de 
Toledo Piza, Phys. Rev. C 53, 334 (1996). 
C.A. Bertulani, V.Yu. Ponomarev, and V.V. Voronov, 
Phys. Lett. B 388, 457 (1996). 

E.G. Lanza, M.V. Andres, F. Catara, Ph. Chomaz, and 

C. Volpe, Nucl. Phys. A 613, 445 (1997). 

P.F. Bortignon, and C.H. Dasso, Phys. Rev. C 56, 574 

(1997). 

V.Yu. Ponomarev, and C.A. Bertulani, Phys. Rev. Lett. 
79, 3853 (1997). 

K. Boretzky et al. (LAND collaboration), Phys. Lett. B 
384, 30 (1996). 

V.G. Soloviev, A. Arima, and N.D. Dang, Phys. Rev. C 
56, R603 (1997). 
V.Yu. Ponomarev, 
Vigezzi, and V.V. 
(1996). 

V.Yu. Ponomarev, P.F. Bortignon, R.A. Broglia, and 
V.V. Voronov, Z. Phys. A 356, 251 (1996). 
C.A. Bertulani, and G. Baur, Phys. Rep. 163, 299 (1988). 
A. Winther and K. Alder, Nucl. Phys. A 319, 518 (1979). 
V.G. Soloviev, Theory of Atomic Nuclei: Quasiparticles 
and Phonons (Institute of Physics, Bristol, 1992). 



P.F. Bortignon, R.A. Brogha, E. 
Voronov, Nucl. Phys. A 599, 341c 



4 



FIG. 1. The energy dependence of 

the ^°®Pb (640-A MeV) + ^°*Pb reaction function calculated 
within first order perturbation theory. The square indicates 
the location of the GDR in ^^^Pb. 

FIG. 2. The cross section for the excitation of the 

2^ component of the DGDR in the reaction ^°*Pb 
(640-A McV) + ■^"**Pb, calculated within second order per- 
turbation theory. The dashed curve shows the contribution 
of the i?l-transition between doorway GDR and DGDR con- 
figurations only. The solid curve is a sum of the above result 
and the contribution of the _El-transitions between complex 
GDR and DGDR configurations. See text for details. 
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